Agent dynamics in kinetic models of wealth exchange 
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We study the dynamics of individual agents in some kinetic models of wealth exchange, particu- 
larly, the models with savings. For the model with uniform savings, agents perform simple random 
walks in the "wealth space". On the other hand, we observe ballistic diffusion in the model with 
distributed savings. There is an associated skewness in the gain-loss distribution which explains the 
steady state behavior in such models. We find that in general an agent gains while interacting with 
an agent with a larger saving propensity. 
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I. INTRODUCTION 

The distribution of wealth among individuals in an 
economy has been a very important area of research in 
economics, for more than a century [H-Q. The same holds 
for income distribution in any society. Detailed analysis 
of the income distribution so far indicate that for 

large income m. 
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where P denotes the number density of people with in- 
come or wealth m. The power law in income and wealth 
distribution is named after Pareto and the exponent ly is 
called the Pareto exponent. The tail of the income dis- 
tribution indeed follows the above mentioned behavior 
and the value of the Pareto exponent ly is generally seen 
to vary between 1 and 3 [1, [Ht3 ■ For any country, it is 
well known that typically less than 10% of the popula- 
tion possesses about 40% of the total wealth and they 
follow the above law, while the rest of the low income 
population, follow a different distribution 3, 5-9]. 

According to physicists, the regular patterns observed 
in the income (and wealth) distribution are indicative 
of a natural law for the statistical properties of a many- 
body dynamical system representing the entire set of eco- 
nomic interactions in a society, analogous to those pre- 
viously derived for gases and liquids. By viewing the 
economy as a thermodynamic system, one can identify 
the income distribution with the distribution of energy 
among the particles in a gas. This has led to several 
new attempts at explaining them, particularly, a class of 
kinetic exchange models have provided a simple mecha- 
nism for understanding the unequal accumulation of as- 
sets. These models are simple from the perspective of 
economics and implement the key factors in socioeco- 
nomic interactions that results in very different societies 
converging to similar forms of unequal distribution of re- 



sources (see Refs. [3,I3, for a collection of large number 
of technical papers in this field) . 

In this paper, we consider the dynamics of individ- 
ual agents. We discuss the money distribution of agents, 
given a particular value of saving propensity. In mod- 
els with distributed savings, we look at a tagged agent, 
and compute the distribution of money gained or lost in 
each interaction. We project the gain-loss behavior into 
a walk in one dimension in a so called "wealth space" . 
Our numerical simulations suggest evidence of ballistic 
diffusion. 



II. GAS-LIKE MODELS 

In analogy to two-particle collision process which re- 
sults in a change in their individual kinetic energy or 
momenta, income exchange models may be defined using 
two-agent interactions: two randomly selected agents ex- 
change money by some pre-defined mechanism. The ex- 
change process does not depend on previous exchanges, 
hence it is a Markov process: 
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where rriiit) is the income of agent i at time t and the 
collision matrix M. defines the exchange mechanism. 

In this class of models, one considers a system with N 
agents (individuals or corporates) and total money M . 
This is a closed economic system where TV and M are 
fixed (microcanonical ensemble), which corresponds to 
no migration or production in the system where the only 
economic activity is confined to trading. Another way 
of looking at this is to consider slow rates of growth or 
decay. Thus the microscopic time scale (of trading) is 
much smaller in comparison to the time scale at which 
the economy experiences growth or collapse. 

In any trading, a pair of traders i and j exchange their 
money |8l-[l3|. locally conserve it, while nobody ends up 
with negative money (rniit) > 0, i.e, debt not allowed): 

mi{t + I) = mi{t) + Am; mj{t + 1) — mj{t) — Am. (3) 
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Time (t) changes by one unit after each trading. 

The simplest case (DY model hereafter) considers a 
random fraction of total money to be shared Q. The 
steady-state {t — ?> oo) money follows a Gibbs distribu- 
tion: P(m) = (l/r)exp(-TO/T); T = M/N. This re- 
sult is robust and is independent of the topology of the 
(undirected) exchange space, be it regular lattice, fractal 
or small- world [llj]. 

Sa ying s is an important ingredient in a trading pro- 
cess [l^l- A saving propensity factor A was introduced 
in the random exchange model [To| , where each trader at 
time t saves a fraction A of its money rriilt) and trades 
randomly with the rest: 

m^{t + 1) = X^n^{t) + e,j [(1 - A)(m,(i) + mj{t))] , (4) 




m,{t + 1) = Am,(t) + (1 - e,,) [(1 - X){m,{t) + m,(t))] , 

(5) 

eij being a random fraction, coming from the stochastic 
nature of the trading. 

In this model (CC model hereafter), the steady state 
distribution P{ni) of money is decaying on both sides 
with the most-probable money per agent shifting away 
from TO = (for A = 0) to M/N as A 1 [l^. This 
model has been argued to resemble a Gamma distribu- 
tion [l3 - [l6| . while the exact form of the distribution for 
this model is still unknown. A very similar model was 
proposed by Angle [l3| several years back in sociology 
journals, the numerical simulations of which fit well to 
Gamma distributions. 

In a real society or economy, saving A is very inho- 
mogeneous. The evolution of money in a corresponding 
model (CCM model hereafter) can be written as 18]: 

'mi{t+l) = \imi{t) + eij [(1 - \i)mi{t) + (1 - \j)mj{t)] , 

(6) 

mj{t+l) = Xjmj{t)+{l-eij) [(1 - Xi)nii{t) + (1 - Xj)mj{t) 

(7) 

It looks similar to the CC model, except that A^ and 
Xj, the saving propensities of agents i and j, are dif- 
ferent. The agents have saving propensities, distributed 
randomly and independently as A(A), such that A(A) is 
non- vanishing as A — > 1; A^ is quenched for each agent (Ai 
are independent of trading or t). The actual asset dis- 
tribution P{m) in such a model depends on the form of 
A(A), but for all of them the asymptotic form of the dis- 
tribution will become Pareto-like (Eqn. (HJ). For uniform 
distribution, A(A) = 1, v = 1. However, for distributions 

A(A)cx(l-A)^ (8) 

P(to) - TO-(2+'5) [lll,[ll[l3. In the CCM model, agents 
with higher saving propensity tend to hold higher aver- 
age wealth, which is justified by the fact that the saving 
propensity in the rich population is always high (20| . An- 
alytical understanding of CCM model has been possible 
until now under certain approximations ^21"], and mean- 
field theory [l^ [22i] suggests that the agent with saving 



FIG. 1: Distribution P{mk\Xk) of money for the tagged 
agent k with a particular value of savings Afc in the CCM 
model with uniformly distributed savings (5 = 0). The data 
is shown for a system of A'^ = 256 agents. 

A possesses wealth to(A) in the steady state, where 

to(A) = C/(1 - A), (9) 

with C cx 1/ log( A^) . The precise analytical formula- 
tion of the above models have been recently considered 
with success, and the results have been derived in most 
cases [1^. There have been recent efforts in analyzing 
the CCM model at a microscopic scale Variations 
of these models have been considered in several manners 
to obtain similar (e.g., for annealed case [1^) or differ- 
ent (as on networks) results [1^ [13] ■ It has also been 
shown that the results are invariant even if one considers 
a cononical ensemble [28| . Some microeconomic formu- 
lations have also proved to be useful [29j] . 

III. STUDYING DYNAMICS OF AGENTS 

In the DY and CC models, agents are homogeneous. 
DY model is nothing but a special case of the CC model 
where A = 0. In these models, looking at individual 
agents and the whole system are equivalent. On the con- 
trary, the presence of the distributed saving propensity 
(quenched disorder) in CCM model gives it a rich struc- 
ture and calls for a careful look at the local scale (at the 
level of individuals) besides computing global quantities. 

In this study, we perform extensive numerical simula- 
tions with system of N agents, with uniform distribution 
A(A) = 1, bounded above by 1 - 1/iV. We look at the 
dynamics of a tagged agent fc, having a saving propen- 
sity Afc, in a pool of N agents distributed according to 
a quenched A(A). We try to see how the individual dis- 
tributions P(7m| Afc) look like (see Fig. [1]). As reported 
elsewhere [lil, [18, 30, 31], the agents with smaller values 
of savings Afc barely have money of the order of average 
money in the market. On the other hand, agents with 
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FIG. 2: Distribution Z)(Amfc|Afc) of money difference Arrife 
for the tagged agent k with a particular value of savings in 
the CCM model with uniformly distributed savings (5 = 0). 
The data is shown for a system of A'' = 256 agents. The inset 
shows that for higher A, probability of losses become larger in 
the region —0.2 < Aruk < 0. 

high saving propensity A^ possess money comparable to 
the average money in the market, and in fact, for the 
richest agent, the distribution extends almost upto the 
total money M. 

A. Distribution of change in wealth 

Upon trading with another agent I, the money of the 
tagged agent k changes by an amount 

Awfe = Tjikit + I) - mk{t) = - (m;(t + 1) - mi{t)) . 

We compute the distributions D{/S.mk\\k) in the steady 
state, given that agent k has a saving propensity Aj, 
(Fig. [2|) . We observe this distribution to have asymme- 
tries both for small and large values of saving propensities 

Total money remains constant in the steady state for 
any agent. An agent with a relatively higher A incurs 
losses which are considerably small in magnitude. This 
immediately suggests that agents with larger savings 
must be having more exchanges where losses, however 
small, occur. The magnified portion of the distribution 
shows that it is really so (shown in the inset of Fig. [2]) . 

B. Walk in the wealth space: Definition 

To investigate the dynamics at the microscopic level, 
one can conceive of a walk for the agents in the so called 
"wealth space" , in which each agent walks a step forward 
when she gains and one step backwards if she incurs a 
loss. The walks are correlated in the sense that when two 
agents interact, if one takes a step forward, the other has 



to move backward. On the other hand, two agents can 
interact irrespective of their positions in the wealth space 
unlike Brownian particles. 

Once the system is in the steady state, we define x{t-\- 
1) = x{t) + 1 if our tagged agent gains money, and x{t + 
1) = x{t) — 1 if she loses. In other words, x{t) performs 
a walk in one dimension. Without loss of generality we 
start from origin (a;(0) = 0), and we insist i = is well 
within the steady state. We investigate the properties 
of this walk by computing the mean displacement {x{t)), 
and the mean square displacement {x'^{t)) ~ {xit))"^ . 

Actually, one can also consider a walk for a tagged 
agent where the increments (i.e., step lengths) are the 
money gained or lost at each step, but the exponential 
distribution obtained for such step lengths (Fig. [2]) indi- 
cates that it will be simple diffusion like 3^]. 



C. Results for the walk 

For the CC model, for any value of the fixed saving 
propensity A, we obtain a conventional random walk in 
the sense {x{t)) is zero and (x2(t))-(x(i))2 - t. However, 
for the CCM model, results are quite different. It is found 
that {x{t)) has a drift, {x{t)) ~ a{Xk)t. a(Afe) varies 
continuously with Afe, taking positive to negative values 
as one goes from low to high values of savings A^ (see 
inset of Fig. [3]) respectively. It is obvious that for some 
Xl, there is no drift, a(A^) = 0. A^ is estimated to be 
about 0.469 by interpolation method. On the other hand 
(x^) — (a;)^ ^ for all A^, which is a case of ballistic 
diffusion (Fig.|l]). 

The negative or positive drifts of the walks indicate 
that the probabilities of gain and loss are not equal for 
any agent in general. Plotting the fraction of times the 
tagged agent gains/loses in Fig. [SJa), it is indeed found 
that an agent with a smaller A gains with more probabil- 
ity while the opposite happens for agents with larger A. 
Indeed, the intersection of the two curves is the point A^. 
where the probabilities are equal and the corresponding 
walk should show (x) = 0. It is difficult, however, to de- 
tect numerically A^ exactly, which lies close to 0.47, and 
check whether an agent with A^ behaves like a conven- 
tional random walker or shows ballistic diffusion. Simu- 
lations using values of A even very close to 0.47 always 
show ballistic behavior. 

In order to explain the above results, we investigate 
at a finer level the walk when the tagged agent with Afc 
interacts with another agent with saving A. First, we 
calculate the average (A) when the tagged agents loses 
or gains and find that for a gain, one has to interact 
with an agent with a higher A in general. This is shown 
in Fig. [5fb). In fact, the average value is very weakly 
dependent on Afc and significantly greater/less than 0.5 
for a gain/loss. This is contrary to the expectation that 
gain/loss does not depend on the saving propensities of 
the interacting agents. 

Having obtained evidence that gain/loss depends on 
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FIG. 3: Measures for the gain-loss walk: the inset shows 
(x) with time for different values of savings A^, showing the 
drifts. The slopes a(Afc) are also shown. The estimate of 
is approximately 0.469. The data is shown for a system of 
N = 256 agents. 



FIG. 5: (a) Plot of probabilities of gain and loss for different 
values of savings propensity A^ of the tagged agent. The data 
is shown for a system of A*" = 256 agents, (b) The plot of the 
average value (A) when a gain/loss is being incurred shown 
against Xk of the tagged agent. 
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FIG. 4: Measures for the gain-loss walk: time variation of 
(x^) — {x)"^ for different values of savings A^, and a guide to 
. The data is shown for a system of A'' = 256 agents. 
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FIG. 6: Data collapse for Xk — 0.2, 0.4, 0.6 and 0.8 is shown for 
Pg — Pi versus the scaled variable y = ^^-^^^^^^ for A'^ = 256. 
The inset shows the unsealed data. 



the interacting agents' saving propensities, we compute 
the probability of gain and loss, Pg and Pi respectively, 
as a function of A for the agent with saving Afc . The data 
shows that indeed an agent gains with higher probability 
while interacting with an agent with A > Afc and vice 
versa. In fact, the data for different A^ collapse when 
Pg — Pi are plotted against a scaled variable y = ^ 
as shown in Fig. [B] indicating a linear variation witfi y, 
i.e.. 

We have checked that there is hardly any finite size ef- 
fect on the collapse in the sense that the similarly scaled 
data for = 100 collapse exactly on those for A^ = 256. 
An agent with a high value of A will interact with a 



higher probability with agents whose saving propensi- 
ties are lesser, causing a loss of money. Therefore in the 
wealth space it will have tendency to take more steps in 
the negative direction. This explains the negative drift 
for large A. 

It is possible to estimate the value of A^ using Eq. ((TU]) 
utilizing the fact that the integrated value of Pg — Pi over 
all A should be zero for A^ = A^. This gives 

1 - (1.5 + 2X1) [log(2.5 + A^.) - log(1.5 + A^)] = 0, (11) 

solving which we get A^ ~ 0.4658 which is consistent with 
the earlier observations. 

It may be added here that in principle the probability 
of gain or of loss while two agents interact, can be calcu- 
lated from the money distribution. In the CCM model, 
when two agents with money mi and m2 and saving 
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propensities Ai and A2 respectively, interact, the differ- 
ence in money before and after transcation for, say, the 
second agent is given by [(1 — Xi)mi — (1 — A2)m2] /2. 
Therefore for the second agent to lose, m2 must be 
greater than m' = "^^i^i^^^ ^ and the corresponding prob- 
ability is given by 

pM r-M 

I P{mi\Xi)dmi I P{m2\X2)dm2. (12) 

Jo Jm' 

However, the exact form of the money distribution is not 
known [2g| for the CCM case. For the CC model, A = 
Ai = A2 and letting M 00, the above integral becomes 

pOO POO 

/ P{mi\X)dTni / P{m2\X)dni2 
Jo Jmi 

POO 

= / P(toi|A)[1 -P(mi|A)]dmi, (13) 

where P{m) = JJ" P{m)dm is the cumulative distribu- 
tion of money. Since F = the R.H.S of Eqn. ^ is 
equal to 1/2 independent of the form of P{m\X). Thus 
in the CC case we find equal probability of gain or loss 
leading to a simple random walk. In the CCM, however, 
the results are expected to be dependent on Ai, A2 as 
indicates. 



IV. DISCUSSIONS 

The analogy with a gas like many-body system has 
led to the formulation of the kinetic exchange models 
of markets. The random scattering-like dynamics of 
money (and wealth) in a closed trading market, in anal- 
ogy with energy conserved exchange models, reveals in- 
teresting features. Self-organization is a key emerging 
feature of these simple models when saving factors are 
introduced. These models produce asset distributions 
resembling those observed in reality, and are quite well 
studied now [llj]. These have prospective applications in 
other spheres of social science, as in application in policy 
making and taxation, and also physical sciences, possibly 
in designing desired energy spectrum for different types 
of chemical reactions. 



In this paper, we have looked at the dynamics of agents 
at the scale of individuals. We study the distribution 
of money for a tagged agent given a particular value of 
saving propensity. We also analyze the distribution of 
money differences in successive exchanges. We conceived 
of a walk in an abstract space performed by the agents in 
different kinetic gas like models which reveals the charac- 
teristics of gains and losses made by the agents. Specifi- 
cally, considering only whether an agent gains or loses, a 
walk can be defined in the wealth space, which is a ran- 
dom walk for the CC model while for the CCM model it 
is found to be ballistic in nature for generic values of A. 

On studying the dynamics at a microscopic level, we 
find that an agent gains with a higher probability when 
interacting with another agent with a larger A. Thus one 
would expect, that for A = 0.5, there will be equal gains 
and losses which would give rise zero drift in the corre- 
sponding walk in the wealth space for the CCM model. 
The value of A^ for which we do get such a result is 
quite close to this estimate. An accurate estimate of 
A^ ~ 0.4658 is obtained by considering the scaling be- 
havior of (probabilities of) gains and losses. 

In the CCM model, our study leads to the discovery 
of the way gains and losses are dependent on the saving 
propensities of the agents, which cannot be arrived at 
using any existing results, either analytical or numerical. 
Using Eqn. (|9]), one can calculate the average money ex- 
changed between two agents which turns out to be inde- 
pendent of their saving propensities. Hence the intriguing 
question that remains to be solved is why the probability 
of gain over loss depends on the savings of the interacting 
agents. 
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